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ABSTRACT 

We present a new general scheme for calculating the structure and dynamics of radiation pres- 
sure driven, photoionized flows. The new method goes one step beyond the Sobolev approxi- 
mation. It involves a numerical solution of the radiative transfer in absorption lines including 
the effects of differential expansion and line interactions such as line locking and blanketing. 
We also present a new scheme for calculating the radiation pressure due to trapped line pho- 
tons in finite, differentially expanding flows. We compare our results for the radiation pressure 
force with those obtained using the Sobolev approximation and show the limitations of the lat- 
ter In particular, we demonstrate that the Sobolev method gives a poor approximation near 
discontinuity surfaces and its neglect of line blanketing can lead to erroneous results in high 
velocity flows. We combine the newly calculated radiation pressure force with self-consistent 
photoionization and thermal calculations to study the dynamics and spectral features of broad 
absorption line flows and highly ionized gas flows in AGN. A comparison with Sobolev-type 
calculations shows that the latter over estimates the flow's terminal velocity and, conversely, 
under estimates its opacity. We also show that line locking on broad emission lines can have 
a significant effect on the dynamics and spectral features of AGN flows. 

Key words: ISM: jets and outflows — galaxies; active — galaxies: nuclei — quasars: ab- 
sorption lines — radiative transfer — methods: numerical 
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1 INTRODUCTION 

Spectroscopic observations of type-I (broad line) active galactic 
nuclei (AGN) show a myriad of UV and X-ray absorption lines 
which, in many objects, are blueshifted relative to the emission 
lines, implying outflows with a wide range of velocities. Such out- 
flows can be crudely divided into narrow absorption line (NAL) 
flows and broad absorption line (BAL) flows. NAL flows are com- 
mon in Seyfert 1 galaxies and are characterized by low velocities 
{<, 1000 kms-i) and narrow absorption profiles (full width at 
half maximum, FWHM~ 100 km s"^). BAL flows are observed in 
~ 10% of luminous AGN (BALQSOs, which also show NAL fea- 
tures) showing sub-relativistic velocities (<J 20, 000 kms~^) and 
FWHM of the same order. The location and mass of the gas respon- 
sible to the observed NAL and BAL flows are poorly constrained 
by observations (e.g., Kraemer et al. 2001, Netzer et al. 2002, de 
Kool et al. 2002 and Everett, Koenigl, & Arav 2002). Theoretical 
modeling is therefore important for understanding the origin and 
the mass flux associated with such flows. 

The modelling of AGN outflows stems from the early studies 
of stellar flows involving acceleration by radiation pressure force 
(Castor Abbott, & Klein 1975, hereafter CAK; but see e.g.. Bland- 
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ford 1994 and for additional types of models). Stellar winds are 
assumed to be in LTE and are driven by radiation pressure in some 
lO"" lines. Most works in this area used the Sobolev approximation 
to calculate the radiation pressure force in the lines (CAK). The 
method has also been used to model AGN flows (e.g., Arav, Li, & 
Begelman 1994, Murray et al. 1995, and Proga, Stone, & Kallman 
2001) despite the very different physical conditions. However, there 
are several complications associated with AGN flows which render 
Sobolev like methods unwarranted. For example, line overlap (line 
blanketing and locking; e.g., Arav 1996) in high velocity flows can- 
not be modeled by Sobolev type methods. A different complication 
arises in NAL flows where the line width is a few times the ther- 
mal width (e.g., Crenshaw et al. 1999), and cannot be treated by 
the Sobolev approximation. Thus, self-consistent solution of AGN 
flows requires more accurate transfer methods. 

There are several other problems unique to AGN flows that 
must be treated in a more self-consistent way. One example is the 
internal radiation pressure resulting from multiple scattering of res- 
onance line photons. The process has been described and analyzed 
by Elitzur & Ferland (1986; hereafter EF86), yet its application to 
realistic flows requires some modification. 

This work is the first in a series of papers that aim to correctly 
model the structure of radiation pressure driven, photoionized flows 
in AGN. In this paper we present a new scheme we have developed 
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and elaborate on the parts that are significantly different from pre- 
vious results. The new method takes into account line overlap and 
uses a better treatment of the radiation pressure force close to dis- 
continuity surfaces. We also revise and improve the EF86 scheme 
for the calculation of internal radiation pressure. The new scheme 
results in different ionization structure, flow dynamics, and differ- 
ent absorption line profiles. The following papers will apply the 
scheme to several important physical situations, such as shielded 
flows, and will investigate specific issues pertaining to the dynam- 
ics, geometry and thermal stability of AGN flows. 

The paper is organized as follows: In section 2 we present the 
new scheme for calculating AGN flows. In section 3 we demon- 
strate the improved capabilities by comparing the results to those 
obtained by the Sobolev approximation and study the dynamics and 
spectral features in NAL and BAL AGN flows. We summarize our 
results in section 4. 



2 METHOD 

In this section we present the new numerical scheme and emphasize 
the improved results for the radiation pressure, radiation pressure 
force, and the absorption line profiles. 



2.1 The equation of motion 



The equation of motion for radiatively driven flows is 

) -9(r), 



(1) 



dv 1 / dPtot 

V— = - Frad -\ 

ar p \ ar 

where v is the flow velocity, r the radial coordinate, p the flow 
density, Prad the radiation pressure force, Ptot the total (gas and 
radiation) pressure, and g the gravity. We assume that contributions 
due to other forces, e.g., drag force, are negligible. This equation 
is solved in conjunction with a continuity condition p = p{r, v) 
which is defined for the geometry under consideration. The impor- 
tant issues of critical points and flow continuity are addressed at the 
end of this section and in appendix B. 



2.2 The radiation pressure force 

We are interested in flows whose dynamics is dominated by radia- 
tion pressure force. To express this force, we make use of the force 
multiplier formalism (see CAK, Stevens & Kallman 1991, and Ch- 
elouche & Netzer 2001; hereafter CNOl) which assumes that 



Frad = HcM, 



where 



H, 



UeO-TLt. 



(2) 



(3) 



is the Compton factor, Ltot is the bolometric luminosity of the 
source, ctt the Thomson cross section, the free electron density, 
and r the radial coordinate. We note that for relatively neutral gas, 
the definition of He is different and must include Rayleigh scatter- 
ing of low energy photons and Compton scattering of X-ray pho- 
tons on bound electrons. Here we focus on ionized gas for which 
equation |3| holds. We also note that in deeper layers of optically 
thick flows, Ltot is smaller than the bolometric luminosity of the 
source due to self-shielding. 

In equation|2| M is the total force multiplier, 



where M\,h is the bound-bound force multiplier, Mbt is the bound- 
free force multiplier, and Ms is the free-free force multiplier. Mc, 
the Compton force multiplier, equals unity by definition. 

Below we define and discuss each of the terms in equation |4| 
in the context of a general flow characterized by its number density 
profile nnir), velocity profile v{r), and a given thermal and ion- 
ization structure. In all cases we assume a typical type-I AGN con- 
tinuum with a "UV bump" and a flat X-ray continuum (aox = 1.4 
and Qa; = 0.9; Le oc E^" , where is defined over the range 
0.1-lOkeV). Solar metallicity gas is assumed throughout. 



2.2.1 Exact and Sobolev-type calculations of Mhb 

Line absorption is sensitive to differential velocities in excess of the 
thermal velocity. We use an exact method to compute the radiation 
pressure force multiplier , AI^^ (r), for a line I of ion at location 
r, where the velocity is v{r): 



0.015n-^-/iAiLfl, 



/oo 
dE'<t,[E' - Ei{l~v{r)/c)\ 
- OO 



T(E',r) 



(5) 



where n ' is the ion number density, /; the oscillator strength, Uth 
the thermal speed for ion Xi, Lei the monochromatic source lumi- 
nosity at the line energy, Ei, and is the line profile. We assume 
pure Doppler line profiles since damping wings have negligible dy- 
namical effects for the densities associated with AGN flows. The 
optical depth, r, at energy E' — Ei{l — v{r)/c) and position r is 
given by 



r{E',r) 



dr^,^iE',r') + J2 / dr^^(E\r') 



(6) 



where 

dr^^iE',r') = 

0.015n^' (r')^^<j}(E' - Ei,(l ~ v(r') /c))dr' , ^'^^ 

Vth 

and drbf is the bound-free optical depth (assumed to vary little over 
the line profile). As evident from equations |6| and the optical 
depth depends on the ionization structure and the velocity profile 
up to location r. This formalism is adequate for calculating the line 
radiation pressure force including the effect of differential expan- 
sion, continuum optical depth, and line interactions (e.g., line blan- 
keting). The total bound-bound force multiplier is given by. 



Mb 



(8) 



X,i.l 



We refer to this formalism as the exact Afbb method. 

A simpler and widely used method is based on the Sobolev 
approximation. This formalism is outlined in the appendix and as- 
sumes that each section of the flow (one Sobolev length across; 
see appendix) is decoupled from the rest of the flow in terms of its 
line optical depth. A useful definition is the local optical depth for 
Thomson scattering (e.g., Arav, Li, & Begelman 1994), 



(9) 



M = Mbb + Mm + Ma + 



(4) 



where Vs is the sound speed. The Sobolev line force multiplier is 
thus a function of the gas ionization level and t (see appendix). 
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which are both local quantities. This dependence on local proper- 
ties is the great advantage of the Sobolev method over exact radia- 
tion transfer methods since the calculation of the radiation pressure 
force is less time consuming. Below we compare the exact and the 
Sobolev methods for various physical situations. 

Previous works in stellar flows have shown that, for LTE con- 
ditions, some 10^ lines are required to fully account for Mbb (e.g., 
Stevens & Kallman 1991). However, the low density AGN flows are 
very different. Our NLTE calculations show that the population of 
most excited levels is small, and that all ions are effectively in their 
ground level. For this reason, only resonance and inner shell lines 
contribute to Mbb. Our code includes some 2000 resonance and 
inner shell lines compiled from Verner, Verner, & Ferland (1993), 
Behar et al. (2001), and Behar & Netzer (2002). In the NLTE condi- 
tions under study, collisional ionization is not very important, and 
the thermal-ionization structure is determined almost completely 
by the shape of the ionizing continuum and the ionization param- 
eter, Ux- The ionization parameter is defined here as the ratio of 
the photon number density in the energy range 0. 1-lOkeV to the 
hydrogen number density. In this notation, broad line region (BLR) 
clouds are characterized by Ux ^ 10"^ while the highly ionized 
X-ray gas is at Ux — 0.1 (CNOl and references therein). 

Stevens & Kallman (1991) find that the most important lines 
contributing to M are optical and near-UV lines (see their table 
1). The main difference between their work and ours results from 
their LTE assumption where the excited states are populated and 
higher transitions contribute significantly to A/bb- In AGN, the 
largest contribution is due to lines with energies around the UV 
bump (10-50 eV; see table0. This is true for a wide range of Ux- 
Thus, the gas dynamics in AGN is sensitive to the luminosity of the 
UV bump which is not directly observed. Our calculations show 
that only a few lines are needed to account for most of the line ra- 
diation pressure for a wide range of Ux and t. Specifically, 90% of 
Mbb in optically thin AGN gas are due to the ten strongest lines. 
As t increases, more lines are contributing significantly to Mbb 
(e.g., some 100 lines are needed to account for 90% of Mbb for 
Ux = 10'^ and t = 10"*). For extremely large t, Mbb < Mbf . 



Line 


n 


Ml 






^x — J^U 


V^^-^bb — 




Hi (A951) 


8 


0.1 




Hi (A973) 


17 


0.08 




Hi (A1027) 


50 


0.07 




Civ (A1549) 


450 


0.01 






11 — 10~2 

Ux — J^'J 


K^^^hh — 




Hi (A973) 


1.5 


0.16 




Hell (A1027) 


4 


0.13 




Hell (A951) 


0.8 


0.12 




Civ (A1549) 


27 


0.02 






11 — 10~1 

Ux — 


( M, , — 

l.^"bb — 




Hell (A256) 


0.8 


0.18 




Hi (A1216) 


0.8 


0.16 




NeVI (A402) 


2 


0.15 




Civ (A1549) 


0.05 


0.05 






11 — 1 


W^^bb — 


1 \ 
^) 


Sxiv (A427) 


0.9 


0.14 




Sixil (A368) 


0.8 


0.12 




Arxvi (A303.5) 


0.3 


0.07 




Civ (A1549) 












Ux = 10 


(Mbb = 


0.4) 


FeXXIV (A192) 


0.7 


0.14 




FeXXIV (A255) 


0.3 


0.09 




FeXXIII(A133) 


1.0 


0.04 




Civ (A1549) 











Table 1. Contribution of the most important lines to the force multiplier as- 
suming a Sobolev t parameter of 10~* (see text and compare with table 1 
in Stevens & Kallman 1991). Also shown is the bound-bound force multi- 
plier, Mbb, for t = 10"* neglecting line blanketing. Note that in all cases, 
lines whose energies are near the UV bump contribute most to the radiation 
pressure force. The contribution of CiV A1549 to Afbb is shown for all 
cases. 



2.2.2 Continuum absorption 
The bound-free force multiplier is 

Mbf(r) = y — / af{E)LE{E)e-^^^'^UE, (10) 

where (j{E) is the bound-free cross-section from all levels (includ- 
ing inner shells) and T{E,r) is defined by equation |6| As shown 
by Arav, Li, & Begelman (1994), and in CNOl, the contribution 
of Mbf to M cannot be neglected for both optically thick and opti- 
cally thin flows. The reason is (e.g., Netzer 1996) that the integrated 
cross section of a bound-free edge ( J asdE) is comparable to that 
of the strongest resonance transition of the same series. Thus, for 
optically thin gas, for ion Xi, 

M^^ ^ {Ei/E,)-"{-Ke''/mc),fi ^ 2 + a 

M,^' " J^^{E/Eo)-'^OEdE - {Ei/Eor^^^o' > 

where a Ei Eq is the continuum slope between the line energy and 
the edge threshold energy, and a is the continuum slope above the 
edge. For our typical type-I AGN continuum, a good approximation 
in the range 10"'' < Ux < W is Mbf/Mbb ^ 0.04t/^"° *''. 



The contribution of free-free processes can be neglected for 
flow densities < 10*^ cm""* (CNOl) and typical infra-red lumi- 
nosities of type-I AGN. 



2.3 The internal radiation pressure 

The total pressure inside the flow is the sum of the gas pressure, 
Pgas, and the isotropic radiation pressure, Prad, which results from 
the scattering of resonance line photons. As shown by EF86, in 
BLR clouds Prad is of the order of Pgas- The numerical scheme 
provided by EF86 takes into account semi-infinite, uniform slabs. 
According to this scheme, 

Prad _ \ ^ STTU^Vth hu 1 — /?e 

Pga. ^ ^ Sn,c^ ksTx e^^l^BT. _ 1 ' *■ 

where Tx is the excitation temperature which is determined from 
the non-LTE calculation and fi^ the local escape probability. AGN 
flows are stratified, dynamically expanding, and have finite total 
optical depth. Hence, /3e must be modified to account for these con- 
ditions. A good approximation, within the framework of the escape 
probability formalism, is 



4 Down Chelouche and Hagai Netzer 



10 



'10" 



10 




ElitzurS Ferland (1986) 



0.2 0.4 0.6 

normalized column density 



0.8 



Figure 1. The escape probability, /3e, as defined by EF86 (dashed curve), 
compai'ed to the one defined in equation fT3l (sohd lines) for static (Ai) = 0) 
and differentially expanding slabs with a unifomi velocity gradient. In both 
cases Dth = 10 km and the total (static) line optical depth is 100. Note 
the steep increase in /3e towards the boundaries. Also note the dramatic 
effect of differential expansion on /3e (the Ai) = 100 kms^i line). In 
this case, the line optical depth is determined by the velocity gradient rather 
than the column density (see text). 



/3e ~ max( 



where 



1 + Tl 



+ 



(13) 



(14) 



T; (r) denotes the optical depth at location r to the two boundaries 
and Ts is the Sobolev line optical depth (see equation Al). This ex- 
pression is valid provided the gas ionization level changes mono- 
tonically and non-local scattering of emission line photons can be 
neglected. FigureQshows the effect of the finite size of the flow on 
Pe, where close to both boundaries photons escape the gas. Dynam- 
ics has a dramatic effect on Pe, as line optical depths are reduced, 
and photons can escape even from the innermost layers of the flow. 
The deviation from the simplified EF86 method, which assumes 
/3e — (1 + Tj"*")"^, can amount to several orders of magnitude in 
Pe. The next section contains a detailed example of the distribution 
of Prad/fgas for Optically thin and optically thick uniform flows. 



2.4 The numerical scheme 

We have developed a numerical procedure for calculating the radi- 
ation pressure force and the global self-consistent structure and dy- 
namics of AGN flows. For given flow properties, v{r), n{r), and 
Ltot , we calculate the ionization and thermal structure of the flow- 
ing gas under NLTE conditions, using ION2002, the 2002 version 
of the photoionization code ION (Netzer 1996). Given the ioniza- 
tion structure, we construct an optical depth vs. velocity grid for ev- 
ery line. The velocity resolution is 0.3«th of the respective ion. We 
divide the flow into thin layers with velocity steps of 0.3uth of the 
heaviest ion. A finer local grid is used in cases where the gas prop- 
erties change on smaller scales (e.g., ionization fronts). For each 
layer, we calculate the corresponding optical depth grids, taking 
into account the local velocity and the continuum opacity, and ob- 



tain Afbb, Mbf and He- Once these quantities are calculated, Fiad 
is obtained and equation Q can be solved. The code optimizes for 
redundant lines whose contribution is negligible either due to the 
low abundance of the relevant ions or due to negligible flux levels 
at the line energy. Given the radiation pressure force, we can move 
on to calculate the self-consistent flow structure and dynamics. 

Our calculations of the radiation pressure force are consider- 
ably different from those used in previous stellar-type works in the 
sense that they do not depend explicitly on dv/dr. Therefore the 
equation of motion does not have the usual form {r^vdv/dr oc 
(r^vdv/dr)"'; see appendix A) and does not necessarily give the 
standard wind solution t; oc (1 — r/ro)^^^. As shown by Poe, 
Owocki, & Castor (1990) for the case of pure absorption, isother- 
mal flows, non-Sobolev solutions are considerably different from 
Sobolev solutions in terms of their dynamical stability and the lack 
of critical points apart from the sonic point. Our model is different 
in several other ways, e.g., we do not assume a constant thermal 
width for all lines and consider ionization and temperature gra- 
dients across the flow. As we show below, the radiation pressure 
force depends also on the velocity of the flow, and is therefore very 
different from that used by Poe, Owocki, & Castor (1990). These 
differences are added to the inherent differences between AGN and 
stellar flows due to the different geometry. Another major differ- 
ence is that stellar winds start subsonically, while the absorption 
troughs of AGN flows imply flows which are supersonic through- 
out. 

In our model, the flow is characterized by its launching point 
To, the density nH(ro), and the initial velocity dq. Given the 
AGN luminosity, Ltot, and the spectral energy distribution (SED), 
we determine the ionization parameter, Ux{ro) at the base of the 
flow. Gravity is expressed using the Eddington ratio, F, where 
g = He/ pV (see CNOl). The first iteration assumes a constant ve- 
locity, V = va, and p oc r~^. We then calculate the radiation pres- 
sure force, solve equationQand obtain v(r). A new density profile, 
nnir), is obtained from the continuity condition. F^ad is then re- 
calculated and the next iteration starts. The process is repeated until 
1% convergence of the flow density profile is achieved. Typically, 
the code takes several hours to run on a PC linux machine for high 
velocity (~ 10* km s~^) stratified flows. 

The scheme presented in this paper assumes vo = Va, i.e., ini- 
tially transonic flows, and is very different from the method used 
in several previous works where the mass flow rate was determined 
from critical point analyses. In a non-Sobolev calculation, the only 
critical point is the sonic point (e.g., Abbott 1982). In AGN flows, 
the sound speed is determined by photoionization and is not re- 
lated to the Keplerian velocity of the gas. The flow is supersonic 
throughout and the mass flow rate is determined solely by the ini- 
tial conditions. More details are given in appendix B. 



3 RESULTS 

In this section we investigate the differences between the exact 
method and the Sobolev approximation. We first study the radia- 
tion pressure force distribution in uniform density clouds and then 
generalize the results to self-consistent BAL and NAL flows. The 
radiation pressure distribution inside optically thin and thick flows 
is also discussed. We choose cases where hydrogen and helium are 
fully ionized in order to allow a meaningful comparison between 
the Sobolev and the exact schemes. 
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Figure 2. An exact absorption line profile (solid line) and a Sobolev line 
profile (dash-dotted line) for a flow with A?; = 100 km and Njj = 
10^^ cm~^. A broadened Gaussian line profile with FWHM=100 km s^^ 
is shown for comparison (dashed line). The Sobolev method produces a 
wrong profile at both low and high velocities. Gaussian line profiles provide 
an even worse approximation to the true profile. The deviation from a flat 
profile in both the exact and the Sobolev method is due to small changes in 
the ionization structure (see also section 3.4). 

3.1 Exact and Sobolev line profiles 

The fundamental difference between the exact method and the 
Sobolev approximation lies in their definition of the line optical 
depth. This is illustrated in the numerical example shown in fig- 
ure |2| corresponding to a slab of thickness Ai?, a column density, 
Nh = lO'^^ cm~^, and velocity interval = 100 km with 
a constant dv/dr across the slab. The resulting line profiles are 
very different. The Sobolev approximation fails to predict the line 
wings since it assumes a constant optical depth across one Sobolev 
length and thus a local, rectangular line profile. The line wings pre- 
dicted by the exact calculations are reminiscent of the Doppler line 
profile which characterizes the local absorption cross-section in the 
flow. For comparison, we also show a broadened Gaussian line pro- 
file with FWHM=100 km s^^. Such line profiles are sometime as- 
sumed in order to calculate the radiation pressure force. This profile 
is different from both the exact and the Sobolev profiles. The ex- 
act calculation shows that the contribution of the line wings to the 
optical depth is substantial. Such line profiles cannot be adequately 
modeled by the Sobolev approximation whose local line profile is 
infinitesimally narrow. Such inaccuracies can significantly alter the 
radiation pressure force distribution inside a flow, and hence the 
flow dynamics. 

3.2 Exact and Sobolev force multipliers 

Figure|3|shows the radiation pressure force distribution for a cloud 
with Ux = 10~^'^ and A^'^f = 10^^ cm~'^ for various values of the 
velocity gradients Av/AR. For Av — 10 km s^^, the solution is 
almost identical to the static solution (see e.g., CNOl), where AI 
monotonically declines into the gas due to the increasing optical 
depth. For larger Av, AI declines from the illuminated surface and 
then flattens as some Sr < AR. This results from the fact that, up 
to some location, the gas differential velocity from its illuminated 
surface is smaller than (uth) and the gas can be considered static. 







gaussian line profile 


exact line profiles < 
















Vv ^Sobolev 
>\ / line''profiles 






/\\Av=10^ km s"^ 


Av=1 l(m s 



'l0' 10° to' 10'° 10" to'" 



5 r (cm) 

Figure 3. The total force multiplier calculated with the exact formalism 
(solid line), the Sobolev approximation (dot-dashed line), and assuming 
broadened Gaussian lines (dashed line). The Sobolev method is a poor ap- 
proximation for M at the illuminated surface but converges to the exact 
result deeper in the cloud. Broadened Gaussian lines (shown for a micro- 
turbulence velocity of 100 km s~^) provide a poor approximation for M. 

Here (vth) is an average over the most important ions which drive 
the gas, and is usually smaller than the sound speed. At larger 5r, 
the flow velocity exceeds (wth) and the optical depth saturates. For 
a uniform cloud this results in a constant optical depth per unit ve- 
locity, and therefore constant M, which is basically the Sobolev 
AI. In all cases and for all Av, the exact force multiplier exceeds 
the Sobolev force multiplier at the illuminated surface by a large 
factor (figure|3}- This is a general behaviour when sharp discontin- 
uous flow surfaces are present. We shall refer to this as the "skin 
effect". The reason is that the Sobolev approximation is insensitive 
to length scales smaller than the Sobolev length (see appendix). For 
AGN flows, the column density of the skin layer is 

Nh ^ rain (vth) ^,W'^^Ua:^ cm~^. (15) 

The first term depends on (fth) which can be set to ~ 13 km s^^ 
(the sound speed at Te — lO^K) and is valid for a wide range of 
ionization parameters since at much higher temperatures, the highly 
ionized gas is driven by heavy metals whose thermal line width is 
small. The second term in this equation reflects the fact that for 
cm , the contribution of lines to the force multi- 
plier in the subsonic part of the flow can be neglected. To conclude, 
flows with rapidly varying properties over small scales, cannot be 
modeled properly by the Sobolev approximation, and the best ap- 
proximation for AI is obtained by assuming static gas. 

Comparison with the exact and Sobolev methods, we also 
show the effect of broadened Gaussian lines on the distribution of 
AI (figure|2). Gaussian models overpredict AI close to the illumi- 
nated surface and under-predict it for the deeper layers. 

While the Sobolev approximation agrees with the exact cal- 
culations beyond Sv = v{rQ + Sr) — v{ro) ~ («th), it cannot 
account for line interaction, such as line blanketing. The effect of 
line interaction is to reduce the force multiplier, hence the calcu- 
lated i\fb? which does not include this effect, is artificially large. 
Figure|4|shows the effect of line blanketing on Mbb for a constant 
density cloud with A'^^f = lO'^'^ cm~'^ and various values of At;. 
We neglect the skin effect and compare the line force multiplier 
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Figure 4. The relative difference between the line force multiplier com- 
puted with the Sobolev approximation and the exact numerical calculations. 
The calculations assume constant ionization structure across a 10^^ cm~^ 
flow. For the faster moving flows, line troughs overlap, line blanketing be- 
comes important and the exact method deviates more from the Sobolev ap- 
proximation. Line blanketing becomes important with increasing Ati since 
more lines overlap and M^b contributes a larger fraction to M. 
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Figure 5. The effect of self-blanketing on M in non-monotonic flows. No- 
tice the large discrepancy between the Sobolev method and the exact calcu- 
lations due to non-local line interaction for Ar > 5 X 10^'' cm. Note also 
the "skin effect" close to 5r ~ 0, and the small but noticeable effect of line 
blanketing for 3 < <5r < 5 X 10^'' cm. 



towards the non-illuminated part of the cloud. The effect of line 
blanketing is small below Aw ~ 3000 km sT^ but becomes sig- 
nificant and reduces M by a large factor at large velocities over 
almost the entire range of Ux (except for extremely ionized gas 
where Mbb <C !)• This is due to the fact that as the velocity range 
increases, more lines are able to interact. Furthermore, as Au in- 
creases, Mbb contributes more to M, and the effect on the gas dy- 
namics is more pronounced. The effect of line blanketing is non- 
local, and as such, depends on the global ionization structure of the 
flow. For this reason, it must be solved self-consistently for each 
flow model (see section 3.2). Our results are in qualitative agree- 
ment with those of Abbott (1982) who investigated line blanketing 
in stellar atmospheres. 

An important case is self-blanketing which occurs in non- 
monotonic flows. We illustrate this by choosing a generic model 
where a uniform density flow, with a column density of 10^^ cm~^, 
is initially accelerated with a constant velocity gradient 2Av/ AR, 
and then decelerated with the same gradient (dashed line in figure 
13 . The exact calculations clearly show a sharp drop in AI beyond 
the deceleration point since lines trace back their absorption pro- 
files and are exposed to lower flux levels. The Sobolev method, that 
is only defined locally, cannot account for such effects and overesti- 
mates M beyond this point. It also fails to correctly predict Mhh for 
dv/dr = since at this point, the Sobolev optical depth diverges 
and Mbb ~ 0. Also marked in figure |5| are the (small) effects of 
blanketing by other lines, and the skin effect close to the illumi- 
nated surface (Sr — 0). 

The locking of absorption lines onto emission lines is another 
potentially important effect (Korista et al. 1996; Arav 1996). Here 
we consider two kinds of line locking: self line locking where reso- 
nance absorption lines are exposed to an enhanced flux level due to 
emission by the same line (e.g., La — La locking), and line lock- 
ing by different lines. The best known example of the second kind 
is the "La ghost" phenomenon in BALQSOs (e.g., Arav 1996). 

We first consider self line locking. All emission line proper- 
ties are taken from the composite spectrum of Telfer et al. (2002). 
The line fist include Mgll A2798, CiV A1549, La, Nv A1240, 



and Ovi A1035 (the contribution of other lines is expected to be 
small for BLR gas) and a constant full width at half maximum, 
FWHM=4000 km s^^ was assumed. This results in Hi La peak 
flux level being 5 times that of the adjacent continuum. In figure 
|6|we show the ratio of radiation pressure force with and without 
line locking. For low ionization flows (Ux ^ 10"''), the radiation 
pressure force is larger by a factor ^ 2 when line locking is taken 
into account. This is mainly due to the locking of La absorption 
on La emission. For larger t, the contribution of La to the force 
multiplier drops since the line becomes optically thick. For moder- 
ate ionization (10~* < Ux < 10^^), locking by higher ionization 
lines becomes important (e.g., CiV A1549). As Ux increases, the 
effect is smaller but extends to larger t. This results from lines being 
more optically thin for larger t. For highly ionized gas {Ux > 0.1) 
the contribution of optical and near-UV lines to Afbb is negligi- 
ble. Obviously, the above results depend somewhat on the assumed 
emission line properties. 

We have examined the Nv A 1240 — La interaction and its 
effect on the radiation pressure force by computing the ratio of 

the radiation pressure force due to Nv A1240 (A^ai24o) on A/bb 
for various uniform density cloud models. When the flow ve- 
locity exceeds 3000 kms~^ (depending on the FWHM of La), 
N V A1240 is exposed to higher flux levels and M^i240 rises. Once 
V <; 5900 km s^^, A/||;Xo drops due to blanketing with La ab- 
sorption line at low velocities. Our results show that the maximum 
relative contribution of N V A1240 to M is can be approximated by 

M^^Zo ^ „ /^N/H\ / 4000 kms-^ \ / EW(LJ \ 

This implies that, for the assumed SED and gas metallicity, the dy- 
namical effect of N V A1240 locking onto La is small. As the radi- 
ation pressure force is roughly proportional to the gas metallicity, 
the effect is larger for overabundant nitrogen. It is also more sig- 
nificant when the ionizing flux is affected by shielding as will be 
discussed in a forthcoming paper. 
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Figure 6. The effect of self line locking on the total (line and continuum) 
radiation pressure force. The ratio of the radiation pressure force when lock- 
ing is taken into account, to the case where it is neglected is shown as a 
function of Ux and t. Note that line locking is most pronounced for large 
velocity gradients (small t) since lines are optically thin. For moderately 
ionized gas the effect is pronounced for larger t. The radiation pressure 
force in highly ionized gas (Ux > 0.1) is not affected by line locking on 
optical and UV emission lines. 




Figure 7. Internal radiation pressure to gas pressure ratio for a typical 
highly ionized gas cloud with a colum density of 10^^ cm~^. For static 
gas ( Av = 0), this ratio reaches a maximum at the center of the cloud, and 
declines to zero towards the non-illuminated boundary. For a finite Av, a 
constant ratio is maintained across a large section of the cloud, and the peak 
/ Pgas is lower. 



3.3 A comparison of the internal line radiation pressure 

As explained in section 2, the distribution of the internal radiation 
pressure can affect the flow structure. Some examples for static and 
expanding uniform clouds are shown in figureQ Our results differ 
significantly from those of EF86 since in our modified procedure 
(equations 1 121 & ll3t . line photons escape the cloud through both 
(all) boundaries. This suppresses Prad close to the cloud's rims. 
Introducing a constant velocity gradient decreases the line opti- 



Av=0 




Figure 8. Radiation pressure to gas pressure ratio for a gas with Ux 



10"'^ and Nh 



10^ 



The profile is non-monotonic due to the 



varying ionization structure. The left hump is mainly due to the scattering 
of C IV A1549 photons and the right hump due to He II A304 line photons. 
Pjad is suppressed once 5v ^ Vs (note the forking of the dashed curves 
from the solid line for Av = 10^, 10^ km s~^). 



cal depth (figure and, for uniform clouds, maintains a constant 
value of Prad across most of the flow. This is shown in figure^for 
t = CO, t = 10"'', and t = 10~*. Larger velocity gradients result 
in smaller Prad- Thus, differentially expanding clouds are less dis- 
rupted by Prad- For optically thin clouds, the peak Prad/Pgas is at 
the cloud centre. 

The above discussion applied to gas that is optically thin to 
continuum absorption. Optically thick clouds show a more intri- 
cate Prad/Pgas profile, as Can be seen from figure |8| Such pro- 
files have, in general, several maxima (two in this case) due to the 
varying ionization structure. In the case shown, the first hump is 
mainly due to the scattering of C IV A1549 line photons. When self- 
shielding becomes substantial, and the ionization structure changes 
significantly, (Nh <l 10^^ cm~^ in figure|8}, the outward optical 
depth in lines of highly ionized species (C IV A1549 in this exam- 
ple) drops, and line photons do not contribute to Prad any more. 
The second hump is due to scattering in lines of lower ionization 
species (in this case mostly He II A304). Velocity gradients (see the 
Aw = 100, 1000 km curves) suppress Prad once 5v exceeds 
the thermal velocity of the ions whose lines contribute most to Prad ■ 
For smaller Sv, the gas is essentially static as far as optical depth 
effects are concerned and Prad is similar to the static case. 



3.4 Self-consistent flow calculations 

In this section we study the self-consistent structure and dynamics 
of AGN flows using the exact and the Sobolev methods. Compar- 
ing the two, we are able to explore the effect of line locking and 
blanketing. We also discuss the dynamical implications of the skin 
effect. 



3.4.1 A linear approximation 

To assess the dynamical effect of line blanketing and locking we 
write the equation of motion for the supersonic section of the flow 
as 
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dv He , . , „ , 
V — = — (M + 5M) , 
dr p 



(17) 



where 5M is a small perturbation i\5M/M\ < 1) representing the 
effect of line blanketing and locking. As a first order perturbative 
analysis we use the unperturbed velocity profile, v{r), to express 
the ratio of the terminal velocities with and without the perturba- 
tion, 



Voo \VooJ / P 



vdv / dr)dr, 



(18) 



where Vod refers to the unperturbed terminal velocity (the mass 
flow rate is determined by the initial conditions of the flow and not 
from critical point analyses which justifies the use of the present 
perturbation analysis). 

The effect of line locking can be written as 



AA/ = {A~l)M(f>,(v,v,), 



(19) 



where <^e is the emission line profile and the characteristic line 
width (~ 0.5FWHM). We assume a Doppler emission line pro- 
file (j)e oc e~" and A = A{Ux,t) is the factor by which 
the force multiplier increases due to locking and is given in fig- 
ure |6| A Sobolev treatment yields the canonical wind solution: 
vdv/dr =const, i.e., v — Vac{i — ro/r)^''^. Assuming that 
Ux =const across the flow (see e.g., the BAL models of Arav, Li 
,& Begelman 1994) then A =const, and the first order correction 
to the terminal velocity is. 



(20) 



This result is insensitive to the exact emission line profile provided 
that its wings decline as fast or faster than v~^. For an ionized 
flow with Vac ~ 1000 kms~^, Voo increases by ~ 5% due to 
line locking. This approximation is adequate for flows whose line 
optical depth are less that unity. It provides a lower limit when lines 
are optically thick since non-linear effects (e.g., due to the increase 
in M and A) become important. 

The effect of line blanketing is non-local by definition. For 
simplicity we assume it becomes important beyond some velocity 
Vb, and cause a decrease in M by a factor B. In this case. 



AM : 



-eK)(l-B)M, 



where &{vb) is the Heaviside function (0{v > vt) 
Vb) = 0). We thus obtain. 



^ ~ 1 - B 1 - for 5oc > «b 



(21) 

1, e(u < 

(22) 



As expected, line blanketing reduces v^c and is more important for 
high velocity flows. For example, an ionized, high velocity (voc ~ 
10* kms'^) flow for which B ~ 0.2 (see fi gureBl results in a 
terminal velocity which is ~ 10% smaller. 



3.4.2 Numerical calculations 

We now study the effect of line blanketing and locking on more 
realistic flows whose structure, dynamics and spectral features are 
calculated, self-consistently, using the algorithm presented in sec- 
tion 2.4. We adopt the model proposed by Arav, Li, & Begelman 
(1994) where the flow consists of small clouds and occupies only a 
small fraction, e, of the volume. The continuity condition is. 




2 

per V 



const, 



(23) 



Figure 9. Self-consistent velocity profiles for flows which are driven a) by 
continuum and broad emission lines and b) by continuum only. The Sobolev 
solution for the first scenario is shown for comparison (dash-dotted line). 
Note tliat line locking on emission lines is dynamically important and in- 
creases the flow terminal velocity by ~ 20%. The effect of line blanketing 
is also important and tends to flatten the velocity profile at high velocities 
compared to the Sobolev solution. 



where p on and e oc v~^. As explained earlier, we consider 
supersonic flows and iterate until the filling factor profile, e{r, v), 
converges to < 1%. In the present calculations the flow is illumi- 
nated by the unobscured AGN continuum and the effects of exter- 
nal shielding will be discussed in a following paper. For illustrating 
the result we chose Ltot = 10*^ erg s~^, ro — 10^* cm, Ux = 
10~^'^, and e(r-o) — 10"'^'^. We also assume that gravity and pres- 
sure gradient force can be neglected. 

Figure|9|shows a comparison of velocity profiles for flows ex- 
posed to continuum with and without emission lines (see section 
3.2). Ux, Ltot, and ro are the same in both cases. As shown, line 
locking increases the final velocity by ~ 20%, which is in agree- 
ment with equation|20| This results from the fact that, in this model, 
O VI A1035, Lc,, N V A1240, and C IV A1549 all contribute signif- 
icantly to the force multiplier. We note, however, that the relative 
importance of line locking depends on the flux in the UV bump, 
where many lines of highly ionized species are located compared 
to the optical and near UV bands. This bump extends over the wave- 
length range of 100 to 400 A and thus a stronger UV bump results 
in a line locking dynamical effect which is relatively less important. 

For comparison we also show in figure|9|the velocity profile 
obtained with the Sobolev approximation for the case where the 
flow is irradiated by a pure continuum emission and for the same set 
of initial conditions. The Sobolev approximation does not include 
line blanketing and, therefore, results in a larger velocity gradient 
and a higher terminal velocity (by ~ 10% in this case). Our cal- 
culations show that using the Sobolev scheme to model the effect 
of line locking will over estimate the terminal velocity by a fac- 
tor of a few since the absorption lines are exposed to an artificially 
enhanced flux level across the entire flow. 

The different velocity profile obtained when line locking is 
significant, results in a smaller flow column density since e (x 
and Nh = / enudr. The lower column density (by roughly 
50% for the model considered here) and the larger flow velocity 
create a shallower absorption trough compared to the case where 
line locking is discarded (figure fTol . 
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Figure 10. The velocity dependent optical depth of CiV A1549 for flows 
accelerated by AGN continuum with and without broad emission lines. 
Broader profiles with smaller optical depths are obtained when fine lock- 
ing is included (figure|9j. Note that the area under the curves is larger for 
the low velocity flow (thick curve) due to the larger column density. Note 
also that the optical depth is increasing with velocity due to changes in the 
ionization of carbon as a result of self-shielding. 
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Figure 11. Velocity profiles for BAL flows with various Eddington ratios, 
r. The r = 0.004 case produces a marginal outflow. All such profiles are 
very different from Sobolev type solution which are monotonically acceler- 
ating with v(r) oc (1 — ro/r)^/^ (see text). 



The exact scheme allows us to study a new family of continu- 
ous flow models: Non-monotonic flows. Such flows have been stud- 
ied in the framework of the cloud model (e.g., Mathews 1990 and 
CNOl) but were excluded from wind calculations on accounts of 
their inconsistency with the Sobolev assumptions. Non-monotonic 
flows have been observed in low luminosity AGN (Ruiz et al. 
2001) and are also predicted by recent hydrodynamical calcula- 
tions (Proga, Stone, & Kallman 2001). In figure^^we show two 
examples of outflows from an object shining well below the Ed- 



dington limit. In both examples Ux 



10" 



We assume that the flows are initially rotating with a Keplerian 
velocity and conserve angular momentum as they accelerate (i.e., 
g{r) oc r~^(l — ro/r)). We study cases where F = 0.006, and 
r = 0.004 which, for the given initial conditions, imply escape 
velocities of 5000 kms~^ and 7500 kms~^ respectively. The 
r = 0.006 model produces an outflow which decelerates beyond 



10" 



10 



10 

in 

Sio" 



10 



r=0.006 



10 

-12000 



r=0.004 



-10000 -8000 -6000 -4000 
V (km s~^) 



-2000 



Figure 12. T(civ A1549) as a function of velocity for the decelerating flows 
shown in figure fTTI Note the large deviation from the CAK solution that 
predicts a constant optical depth as a function of velocity (see text). 



r ~ 2ro and maintains a coasting velocity of 10, 000 km s~^. The 
r = 0.004 model corresponds to a marginal outflow which barely 
escapes the system (v{r — oo) ~ 0). In such flows, BAL-type ve- 
locities can be obtained for part of the time but the terminal veloc- 
ity can be much lower. Thus, BAL flows can be observed in objects 
which shine at only a small fraction of their Eddington luminosity. 

The velocity dependence of r(ClV A1549) for the above de- 
celerating models are shown in figure fT2l This again is very dif- 
ferent from the results for the monotonic BAL flows considered by 
Arav, Li, & Begelman (1994; see their p oc models), and from 
stellar winds profiles which are characterised by t{v) =const. 

To illustrate the dynamical importance of the skin effect we 
study the dynamics of low velocity flows. We do so by ignoring 
the contribution of pressure gradients to the total force. We assume 
Ltot = 10^^ erg s~\ ro = 10^^ cm and [/^(ro) = 0.1, typical 
of highly ionized gas. As shown in fieure fT3l the exact solution, 
that includes a larger force multiplier close to the illuminated sur- 
face, makes a large difference to the terminal velocity of such flows. 
When the flow is optically thick, the dynamical effect can be more 
dramatic due to the non-linear nature of the opacity. When initially 
subsonic flows are considered (see appendix B for the case where 
p oc r'^v"^), the skin effect will decrease the mass flow rate since 
larger density gradients are needed in order to balance the larger 
radiative force close to the illuminated surface. 

The effect of P^ad on the flow dynamics is small for flows that 
are not very optically thick beyond the He II edge since Piad 
Pgas (EF86). The effect can be neglected once the lines become 
optically thin due to expansion. 



4 CONCLUSIONS AND SUMMARY 

We have presented a new method for calculating BAL and NAL 
flows with different properties. The main improvements over pre- 
vious methods are the non-Sobolev treatment of the lines and the 
treatment of the internal radiation pressure. We show that Sobolev 
based calculations suffer from several limitations and disagree with 
the exact calculations in both high and low velocity flows. The main 
differences are due to line blanketing in high velocity flows and in- 
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Figure 13. Velocity profiles for highly ionized NAL flow {Ux = 0.1) for 
exact (solid line) and Sobolev (dashed-dotted line) calculations neglecting 
pressure gradient effects. Including the skin effect can considerably enhance 
the terminal velocity of such flows. 



correct calculation of the line transfer near discontinuity surfaces 
in low velocity flows. 

We present a new set of models pertaining to non-monotonic 
flows. Such flows cannot be modeled by the Sobolev approxima- 
tion, but are important in astrophysical systems (e.g., Ruiz et al. 
2001). The calculations show that, under some geometries, it is pos- 
sible to obtain high velocity flows in objects that shine at only a few 
per cent of their Eddington luminosity. Such non-monotonic flows 
can reach BAL velocities yet they barely escape the system. 

We also study the effect of radiation pressure due to trapped 
line photons in optically thin and optically thick differentially ex- 
panding media. In such cases, the radiation pressure distribution is 
considerably different from that in static uniform media. Radiation 
pressure is suppressed close to the flow boundaries and can obtain 
several local maxima. Differential expansion tends to suppress and 
homogenize the internal radiation pressure and reduce the dynami- 
cal effects on the flow. 

This work is another step towards constructing more realistic 
dynamical models for radiatively driven flows in AGN. It allows 
modeling of radiative transfer in non- stationary situations and pro- 
vides a mean for incorporating dynamics into realistic photoioniza- 
tion models. The use of detailed radiative transfer codes to model 
such flows is a formidable, yet highly desirable task, as it may break 
the degeneracy between density and location in pure photoioniza- 
tion modelling of AGN where Ux oc Ltotr~^n'^ . This will allow, 
eventually, to constrain the mass entrained in such flows with im- 
plications to the understanding of AGN and their interaction with 
the host galaxy and the inter-galactic medium. 

We thank N. Arav for many enlightening discussions. We ac- 
knowledge financial support by the Israel Science Foundation grant 
no. 545/00. 
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APPENDIX A: THE SOBOLEV FORMALISM 

The Sobolev approximation applies to situations where the flow 
properties change on scales larger than the velocity scale, e.g., 
Vgdr/dv <C pdr/dp, where Vsdr/dv is the "Sobolev length". The 
Sobolev approximation redefines the line optical depth so that each 
section of the flow, within one Sobolev length, is assumed to be in- 
dependent on other sections. The Sobolev line optical depth (which 
replaces equation|6|in this paper) is defined as 



= 1. 



(Al) 



In the Sobolev approximation, the bound-bound force multiplier, 
M^f , is given by. 



Ml 



E 



E,L 



(?) 



(A2) 



This definition of M^^' is somewhat different from that of previ- 
ous works (e.g., Arav, Li, & Begelman 1991, Stevens, & Kallman 
1990; see however Proga, Stone, & Kallman 2000) because we in- 
clude the bound-free optical depth term, Thi{E, r). This is done in 
order to allow a meaningful comparison with the exact numerical 
method (equations|5|&|6j. 

Ignoring the effects of self-shielding by bound-free processes, 
CAK provide a most convenient approximation for M^^, 
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(A3) 



where q ~ 0.9 for AGN (Murray et al. 1995), and fc is a normaliza- 
tion constant. This parameterization was later refined by Owocki, 
Castor, & Rybicki (1988) to prevent the divergence of eguation lABI 
at t <C 1. The CAK parameterization, when plugged into the equa- 
tion of motion, requires vdv / dr =const, which leads to the canon- 
ical stellar wind solution, v = v,^ ( 1 — ro /r) 
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APPENDIX B: CRITICAL POINTS IN A NON-SOBOLEV 
FORMALISM 

Here we follow the arguments outlined in Abbott (1984) and apply 
them to AGN flows. For a spherical continuous flow (p oc r~^t;~^), 
the equation of motion (equation^ can be rearranged to give, 

Idv _ 1 

V dr — Vg 

where niH is the proton mass. This equation has a singular point 
ai V = Vs and a transonic solution can be obtained only if the 
numerator satisfies 

2vl ^ otU.. f _ 1 \ ^ Q ^32) 
r A-Kr^mHC \ r(r) J 

at the sonic point. Equation IB2I has one or more solutions if 
M{r) < l/r(r). This condition is met for line driven winds and 
sources that shine well below their Eddington limit (e.g., stars). 
However, the Eddington ratio in luminous AGN is of order 0.1 
and, unless the flow is highly ionized, A/bf + A^bb > l/F (e.g., 
CNOl and the calculations presented in this paper). Thus, a sta- 
tionary (non-accelerating) hydrostatic subsonic structure cannot be 
formed. Furthermore, observations show that the gas temperature 
in AGN is not related to its location (e.g., the broad line region gas 
whose sound speed is ~ 10 km s^^ moves in Keplerian orbits with 
velocities in excess of 1000 kms"^). 

The above considerations suggest that, in AGN flows, the sub- 
sonic gas is never observed and therefore can be ignored. This 
means that the mass flow rate is determined by the initial condi- 
tions at a location where v > Vs- Physical scenarios which could 
motivate this assumption are the evaporation of cool gas to super- 
sonic velocities from large condensations in the BLR, the central 
accretion disk, the torus, etc. 
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Elitzur & Ferland (1986) 



0.2 0.4 0.6 0.8 1 

normalized column density 
Figure 1: Escape probability for a static medium showing the corrected Elitzur & Ferland approximation 

Erratum: Non-Sobolev modelling of radiation pressure driven flows 
in Active Galactic Nuclei. 

Doron Chelouche and Hagai Netzer 

School of Physics and Astronomy and the Wise Observatory, The Beverly and Raymond Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel 



The above paper (Chelouche & Netzer 2003) includes an erroneous reference to the work of Elitzur & Ferland (1986) which resulted in an 
erroneous Fig. 1. Contrary to our statement in section 2.3, our approach to the calculation of the escape probability in the case of a static 
medium is in fact similar to the one proposed by Ferland & Elitzur (1984) and Elitzur & Ferland (1986). In both methods the escape 
probability, /Jf''"'', increases toward the rims of the cloud (the Elitzur & Ferland expression is = 0.5 [/3f ''"''(r+) + Pf^^'^ir')] 

where are the measured optical depths from points inside the cloud to its rims). The only difference is in the exact way the radiation 
pressure drops toward the edge of the cloud. Both methods are based on the escape probability approximation and we do not consider them 
to be significantly different (this is illustrated in the corrected Fig. 1 above). We apologize for the error and thank G. Ferland and M. Elitzur 
for pointing it out. 
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